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We study the masses and wave functions of heavy-light quark-antiquark system in the relativistic 
potential model, where we focus on the higher excited states of B and D mesons. The resonances 
of Dj (2740), D}(2760), D}(3000) and D*j( 2860) recently observed in experiments are considered. 

We find that the four resonances can be identified as the |1 1,3 D2), 11 3 Z?3), |1 1,3 £3) and 11 3 Z? 3 ) states 
due to the prediction to the states with l = 2 and l = 3 in the potential-model calculation. The 
masses of n = 2 excited states are also calculated. All the observed mesons can be accommodated 
successfully in the potential model. 
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I Introduction 


Heavy-light quark-antiquark system Qq plays an im¬ 
portant role in studying the strong interactions between 
quark and antiquarks. Great progress was made in 
the measurement of the spectrum of heavy-light quark- 
antiquark system in the last decade. In the charm sec¬ 
tor, several new excited meson states were discovered 
in the previous several years [l]-Q . For Dj mesons, 
the resonances of £>(2740)°, £>*(2760)°’+ @, £>*(2650) 
and ££(3000) @| were observed, while for D sJ mesons, 
the resonances £>^(2632) 0, X(2690), £>+ 7 (2860) 0, 
£>jj(2700) 0 and ZJ+j (3040) 3] were established. The 
spectroscopy provides a powerful test of the theoretical 
predictions based on the quark model in the standard 
model. 


The bound state of the heavy-light quark-antiquark 
system can be studied with the relativistic wave equa¬ 
tion, where an effective potential compatible with QCD 
is taken 0, 5]- The potential shows a linear confining 
behavior at long distance and a Coulombic behavior at 
short distance. The properties of B and D mesons have 
been investigated in the relativistic potential model in 
Refs. I8h14||. In our recent work |15|, the spectra for 
lower orbital and radial excitations of B and D mesons 
were calculated, a good agreement with experiment [Iff] 
was obtained. Here we extend our previous work to in¬ 
clude the higher excited states and compare the predic¬ 
tion with the spectrum measured in experiment in recent 
years. The excited states of the heavy-light system have 
also been studied in Ref. 0 in the relativistic quark 
model, which includes the leading order corrections in 
1 /mc,b expansion. The Hamiltonian used in Ref. 17] is 
a kind_of relativistic Dirac Hamiltonian. Compared with 
Ref. 0 , the method used in this work is different and 
the recently observed states are also carefully studied in 
this work. 


The organization of this paper is as follows. In section 
II, we give the effective Hamiltonian for the heavy-light 
quark-antiquark system and solve the relativistic wave 
equation theoretically. Section III is for the numerical 


result and discussion. In section IV we have a brief sum¬ 
mary. 


II The effective Hamiltonian and the solution of the 
wave equation 

In the potential model one assumes that hadrons can 
be treated in terms of valence-quark configurations. Here 
we employ the relativistic potential model to describe 
the heavy-light B and D mesons. The effective potential 
between the quark and antiquark in the meson is one- 
gluon-exchange (OGE) dominant at short distances and 
with a linear confinement at long distances. Since there 
is a light quark in the heavy-light system, it requires to 
include relativistic dynamics in the Hamiltonian. The 
relativistic version of Schrodinger wave equation describ¬ 
ing the heavy-light system is written as 

Hip(r) = Eip(f). (1) 

In the rest frame of the bound state system, the eigen¬ 
values of the Schrodinger equation will be the masses of 
mesons. 

The effective Hamiltonian can be written as 

H = H 0 + H ', (2) 

with 

#o = \Jp{ + ml + \Jj% + m| + V(r), (3) 

where mi, p\ are the mass and momentum of the heavy 
quark, respectively, and m 2 , P 2 the mass and momentum 
of the light antiquark. V(r) is the funnel potential, which 
includes only the spin-independent part of the strong in¬ 
teraction between the heavy and light quarks 

V(r) =--^-L+br + c. (4) 

3 r 

The potential is taken as a combination of a Coulomb 
term and a linear confining term, whose behavior is 
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compatible with QCD at both short- and long-distance 
[ 1 > EM Gil • The first term is obtained from the calculation 
of one-gluon-exchange diagram in perturbative QCD. It 
dominates at short-distance. The second term is the lin¬ 
ear confining term. The third term is a phenomenolog¬ 
ical constant, which can be adjusted to give the correct 
ground state energy-level of the quark-antiquark system. 

The running coupling constant a s (r) in eq.Q is ob¬ 
tained from the coupling constant in momentum space 
a s (Q 2 ) after the Fourier transformation. It can be writ¬ 
ten in a more convenient form Q 



where and 7 \ are free parameters which can be fitted 
to make the behavior of the running coupling constant at 
short distance be consistent with the coupling constant 
in momentum space predicted by QCD. In this work, we 
take ai = 0.15, a 2 = 0.15, <23 = 0.20, and 71 = 1/2, 
72 = vTO/ 2, 73 = v / 1000/2 [3. 

The second term H' in eq.0 is the spin-dependent 
part of the Hamiltonian 

H' = H hyp + H so , ( 6 ) 

where H hyp is the spin-spin hyperfine-interaction term, 
which is 

Q0 7 r 

-ft -aj C) 

3 777-1777-26 V 6 \ J 

with 

Sa(r) = (^) 3 e- CTV , (8) 

y/TT 

where the parameter a is taken as quark mass-dependent 

i 


<7 = 



1 / 4toiTO 2 \ 4 \ 2 / 2mim 2 \ 2 

2 V( TO i + m 2) 2 ) J +S ° + m 2 ) 

(9) 


here <jq and so are phenomenological parameters. 

The spin-orbit interaction term H so consists of the 
color-magnetic term and the Thomas-precession term 


rso 


4 a sir ) / J_ + _J_\ ( si • L 
3 r 3 \mi fri 2 cj l JTii 



1 dV(r) { si ■ L s 2 -l \ 
2 r dr y m 2 (m 2d ) 2 J 


( 10 ) 


The spin-dependent Hamiltonian H' is the result of 
OGE diagram in the non-relativistic approximation [3,0] ■ 
It is reasonable that there might be contributions of non- 
perturbative dynamics in the bound state system and 
relativistic corrections for the light quark. In this work, 
the S(r) in the spin-spin contact hyperfine interaction 
term is replaced by (-2 ='j 3 e ~ <j2r2 as in Ref. [20]. The 
mass of the light quark 7712 ’s in the denominators of eqs. 
0 and m are replaced by a set of new parameters 
fri 2 i , i = a, b , c, d, which is assumned to include the 
relativistic corrections and the bound-state effect in the 
heavy meson [l5|. In the original expression of the spin- 
dependent interaction, the momenta of the quarks are 
dropped. However, it is not reasonable to drop them 
completely for the heavy-light system as the light quark 
should be highly relativistic. There may also be bound- 
state effect in the spin-dependent interaction terms which 
can not be described by the OGE potential. We assume 
that these effects can be included by replacing the light 
quark masses with a set of new parameters, which can 
be determined by fitting the spectra of heavy mesons. 
One can find in the next section that this assumption 
does work, all the masses measured in experiment can be 
accommodated well. 

In our previous work 15(, only the states of B, D, B s , 
D s mesons with lower orbital quantum numbers l = 0,1 
are considered, the masses and wave functions of the cor¬ 
responding states with J p = 0 _ , 1”, 0 + , 1 + , 2 + are 
obtained. In this work we extend our previous work by 
studying more excited states with higher orbital quantum 
numbers. In the heavy quark limit the spin of the heavy 
quark decouples from the interaction 21-29j], therefore 
the spin-dependent interaction can be treated perturba- 
tively. We choose the eigenstates of spin-independent 
Hamiltonian Hq as the basis of perturbative expansion. 
The eigenstates of H 0 is denoted as \JM,sl), where J 
is the quantum number of the total angular momentum, 
M the magnetic quantum number, s the total spin, l 
the relevant orbital angular momentum. The 5 a (r) term 
of the hyperfine interaction H hyp conserves all the four 
quantum numbers, it gives the splitting of states with 
different total spin. The tensor part of H hyp does not 
conserve the orbital angular momentum, it causes mix¬ 
ing between the states with different orbital angular mo¬ 
menta 3 Lj - 0 - 3 L'j, while the spin-orbit interaction H so 
does not conserve the total spin, it causes mixing be¬ 
tween the states with different total spin quantum num¬ 
bers 1 Lj -O 3 Lj. 

The details of solving the wave equation have been 
~|, which is omitted in the present paper 


14, 


given m 

for brevity. The readers can refer to these references. 

The mass matrix is calculated perturbatively in the 
basis of \JM, si ), which can be written as 


where L = rxP is the relative orbital angular momentum 
between the quark and antiquark. 


(Hu H12 \ 

V H 2 i H22 J ’ 


( 11 ) 
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where ffy = (ipi\H\ipj), with i,j = 1, 2, tpij’s denote the 
eigenstates of the Hamiltonian Hq . By diagonalizing the 
mass matrix, we can get the masses and wave functions 
of the mixing eigenstates. The perturbative contribution 
of the spin-dependent Hamiltonian H' to the eigenvalues 
of states with high orbital quantum numbers are given 
below (Contributions to states with lower orbital mo¬ 
mentum have been given in our earlier work, Ref. 0). 

(1) The mass matrix of state J p = 2~ 

Both states with s = 0, l = 2 and s = 1, l = 2 can 
construct the J p = 2~ state. The basis for the mixing is 
\ipi) = 1 1 Z? 2 ) and ^ 2 } = l 3 ^)- The matrix elements of 
the mass matrix are 

Hu = ( 12 ) 

H i2 = \l^(^i= 2 {r)\h 1 (r)-h 2 {r)\^%(r)), (13) 

H 21 = H* 12 , (14) 

H 22 = + + 

~\hi{r) ~ \h 2 {r)ip\%{r))). (15) 

With eqs. 0 ~ (HU), the cases with more j 1 !^} and 

\ 3 D 2 ) mixing states can be obtained. 

(2) The mass matrix of state J p = 3“ 

The J p = 3“ state is mixture of 3 D 3 and 3 G 3 states, 
both states with s = 1, l = 2 and s = 1, l = 4 can 
construct the J p = 3“ state. The basis for the mixing is 
= \ 3 D 3 ) and 1 ^ 2 } = | 3 Gs)• The matrix elements of 
the mass matrix are 


Hu = El% + (ri%(r)\[±f(r)-±g(r) 

+hi(r) + h 2 (r)]\^%(r)), (16) 

H 12 = ^\/3(V> ; ( = 2 (r-)|5WlW=4( r ))> ( 17 ) 

H 2 1 = H* 2 , (18) 

H 22 = E[% + (if>i= 4 (r)\[^f(r) - ^g(r) 

- \h 2 {r)]\tp\%(r)). (19) 


With eqs. m~m, the cases with more | 3 Z? 3 ) and 
| 3 G 3 ) mixing states can be obtained. 

(3) The mass matrix of state J p = 3 + 

Both states with s = 0, l = 3 and s = 1, l = 3 can 
construct the J p = 3 + state. The basis for the mixing is 
|'i/q} = l 1 ^) and \ij} 2 ) = l 3 !^). The matrix elements of 


the mass matrix are 

Hu = -Bj= 3 - |(V’j ( =3( r )l/( r )lW=3( r ))> ( 2 °) 

H 12 = - h 2 (r)\i/jl%(r)), ( 21 ) 

H 21 = H* 2 , (22) 

H 22 = E\% + {^\%(r)\^f(r) + ^g{r) 

-^iW -^2 (r)i>i%(r))). (23) 

With eas. (EUl) ~ (B^l) . the cases with more | 1 Z? 2 ) and 
1 3 D 2 ) mixing states can be obtained. 

(4) The mass matrix of state J p = 4 + 

The J p = 4 + state is mixture of 3 if 4 and 3 H 4 states, 
both states with s = 1, l = 3 and s = 1, l = 5 can 
construct the J p = 4 + state. The basis for the mixing is 
\ipi) = | 3 ^ 4 ) and |^ 2 ) = l 3 ^)- The matrix elements of 
the mass matrix are 


H u 

= E^% + {ipt%(r)\[-f(r) - -g(r) 



+ (r) + ^h 2 (r)]\ipl%(r)), 

(24) 

h 12 

= -y W’l= 3 {r)\g(r)\^%(r)), 

(25) 

H 2 i 

_ rr* 

(26) 

h 22 

= El% + (4%(r)\[\f(r)-±g(r) 



—3hi(r) - Sh 2 {r)]\ijjl%(r)). 

(27) 


With eqs. (HU) ~ J23, the cases with more j 3 ^) and 
| 3 ii 4 ) mixing states can be obtained. 

(5) The mass matrix of state J p = 4~ 

Both states with s = 0, l = 4 and s = 1, l = 4 can 
construct the J p = 4~ state. The basis for the mixing is 
|i pi) = | 1 G 4 ) and \i/j 2 ) = | 3 G 4 ). The matrix elements of 
the mass matrix are 

Hu = ^i =4 - ^(V'*= 4 (»")I/(»")IW= 4 (’")) > ( 28 ) 

H 12 = V5(^i=4(r)\hi(r)-h 2 {r)\i>l%(r)), (29) 

H 21 = H* 2 , (30) 

H 22 = E[% + {^ 4 (r)\^f(r) + ^g(r) 

~\h! (r) - ^h 2 {r)i>[%{r))). (31) 

With eqs. ED ~ (ED, the cases with more | 1 G 4 ) and 
| 3 G 4 ) mixing states can be obtained. 

( 6 ) The mass matrix of state J p = 5 + 

Both states with s = 0, l = 5 and s = 1, l = 5 can 
construct the J p = 5 + state. The basis for the mixing is 
\ipi) = I 1 # 5 ) and \ijj 2 ) = | 3 ^ 5 )- The matrix elements of 
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the mass matrix are 

Hn = Ej%-^%(r)\f(r)\ri%(r)), (32) 

h 12 = \j^(^%(r)\h 1 (r)-h 2 {r)\^%(r)}, (33) 

Hi i = HI 2 , (34) 

iJ 22 = E[°{ + <^ ( = 5 ( r )l^/( r ) + 7}ff( r ) 

~\hi(r) - \h 2 {r)if}[%(r))). (35) 

With eas. (l32l) ~ (l35l) . the cases with more and 

\ 3 Hs) mixing states can be obtained. 

(7) The mass matrix of state J p = 5~ 

The J p = 5“ state is mixture of 3 Gs and 3 /5 states, 
both states with s = 1, l = 4 and s = 1, l = 6 can 
construct the J p = 5~ state. The basis for the mixing 
is |V>i) = | 3 G 5 ) and |^> 2 ) = | 3 T 5 ). The matrix elements of 
the mass matrix are 


H u 

= e\% + (vSMIlJ/M - ^g(r) 



+2ht(r) + 2 /i 2 (r)]|V>}° 4 (r)), 

(36) 

h 12 

= YY\/y' (V’ElWbWIV’EeW). 

(37) 

h 21 

_ /T * 

^125 

(38) 

h 22 

= £/ ( °6 + (WieMItj/W - ^jff(r) 



~\ h !(r) - \h 2 {r)}\il)\%{r)). 

(39) 


With eqs.® ~ flMl) , the cases with more | 3 G 5 ) and | 3 / 5 ) 
mixing states can be obtained. 

In the above equations, the functions f(r), g(r), h\(r) 
and h 2 (r) are defined as 


f(r) 

9(r) 

hi(r) 

h 2 [r) 


o--— a s (r)6 a (r), 

9mim 2a 


4 a s (r ) 1 
3 m\fh. 2 b r 3 ’ 


"4 a s {r) 

( 1 

+ 1 ^ 

3 r 3 


m 2c J 


1 dV{r) 1 ' 

1 


2 r dr m\ 

m i ’ 

’4 a s (r) / 1 

h M 1 

3 r 3 \mi 

rh 2c J rh 2c 

1 dV(r) 1 




2 r dr ( m 2d ) 2 


(40) 

(41) 


(42) 


(43) 


III Numerical result and discussion 

The parameters used in this work are the quark masses, 
the potential parameters b,c,rh 2 i, cro and Sq. Compar¬ 
ing with our previous work [151 ] . we have slightly adjusted 


some of the parameters. The ground states are not sen¬ 
sitive to the value of the parameter 6 , which is slightly 
increased to give the correct predictions of excited states. 
The values we obtain by fitting are 


m b = 4.99 GeV, m c = 1.58 GeV, 
m s = 0.32 GeV, m u = rrid = 0.06 GeV, 
b = 0.175 GeV 2 , c =-0.312 GeV, 
cTo = 1.80 GeV, So = 1.55. (44) 

The values of rh 2a , rh 2 b 1 rh 2c and fh 2 d depend on the 
quark-antiquark system, they can be written as 


m 2 i — eifh 2 , i = a,b, c, d. 

We find the values of e*’s and m 2 are 

(e a ,e b ,e c ,e d ) = (1.10,1.30,1.30,1.32) 
for (bq) and icq) systems, and 

(ea,eb,e c ,e d ) = (1.10,1.10,1.10,1.31) 
for (bs) and (cs) systems, and 


(45) 


(46) 


(47) 


m 2 = 


0.562 GeV for {bq) system, 
0.679 GeV for ( 6 s) system, 
0.412 GeV for (cq) system, 
0.488 GeV for (cs) system, 


(48) 


here q is the light quark u or d. 

Numerical calculation shows that the solution is stable 
when L > 5 fm, N > 50. Here we take L = 10 fm, 
N = 100. The definition of L and N can be found in 
Refs. [H,[l5|. 

The radial wave functions 'tp n i{r) for D meson in co¬ 
ordinate space are depicted in Fig.l as an example. 
\4’ni{'r)\ 2 r 2 is the possibility density distributed along 
the quark-antiquark distance r. We can use the Fourier 
transformation of the wave function in the coordinate 
space to get the momentum distribution, which is useful 
for studying B and D decays. 

In this work we calculate B, D , B s , D s bound states 
with higher orbital quantum numbers, the quantum num¬ 
ber J p of the new corresponding states are 2“, 3 _ , 3 + , 
4 + , 4~. the numerical results for B , D , B s , D s bound 
states are presented up to l = 4 in four tables. The nu¬ 
merical results for D, D s bound states are given in Table 
I and II, the numerical results for B , B s bound states 
are given in Table III and IV. For completeness, results 
for ground-states and lower orbital excitations are also 
given. The mass predictions for lower orbital excitations 
have been calculated in Ref. [ 8 ] several decades before, 
which are also presented in Tables I to IV for compari¬ 
son. The masses for most of the lower excitation states 
obtained in this work are approximately consistent with 
the masses given in Ref. jsj]. Only for a few mesons, such 
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FIG. 1: The radial wave function i/jni(r) of D meson in co¬ 
ordinate space, l from 0 to 4. (a) is for n = 1, (b) is for 
n = 2. 


as £i(2420), £*s0(2317), £*sl(2460) etc., the mass pre¬ 
dictions are improved, the results in this work are more 
consistent with experiment. Mixings between states with 
appropriate quantum numbers are considered in our cal¬ 
culation. For doublets with J p = 1 + , 2 _ , 3 + , 4~, mixing 
occurs more strongly. The mixing from state of higher 
radial quantum number with n = 2 can not be neglected, 
it can be as large as 10%. The masses measured in exper¬ 
iments are given in the last column in Tables I and III, 
from which we can see that the theoretical calculation in 
this work can accommodate the experimental data well. 


Many new resonances with open charm flavor were ob¬ 
served in experiments recently Si- Our mass predic¬ 
tions for these resonances are consistent with the experi¬ 
mental data. The masses measured by the LHCb collab¬ 


oration are listed in the following 


M(£>i(2420)°) 
M(D%( 2460)°) 
M(D}( 2650)°) 
M(D}( 2760)°) 
M(£j(2580)°) 
M(Dj( 2740)°) 
M(£j(3000)°) 
M(D* 2 { 2460)°) 
M(£}(2760)°) 
M(£}(3000)°) 
M{D* 2 (2460)+) 
M(£}(2760)+) 
M(£}(3000)+) 


2419.6 ± 0.1 ± 0.7 MeV 

2460.4 ±0.4 ±1.2 MeV 
2649.2 ±3.5 ±3.5 MeV 

2761.1 ±5.1 ±6.5 MeV 

2579.5 ±3.4 ±5.5 MeV 
2737.0 ±3.5 ±11.2 MeV 
2971.8 ±8.7 MeV 
2460.4 ±0.1 ±0.1 MeV 

2760.1 ±1.1 ±0.1 MeV 

3008.1 ±4.0 MeV 

2463.1 ±0.2 ±0.6 MeV 

2771.7 ±1.7 ±3.8 MeV 

3008.1 MeV. 


The assignments of the above observed states are listed 
in Table HI and ITT1 Dj( 2740) can be identified as mixing 
state of ll 1 ^) and |1 3 1?2) with J p = 2~. £}(2760) can 
be identified as dominantly |1 3 £ 3 ) state with J p = 3“, 
which agrees with Ref. [33 ■ £}(3000) can be identified 
as mixing state of and |1 3 ±3) with J p = 3 + . The 

other three observed states are £j(2580), £}(2650) and 
£j(3000), are identified as n = 2 excited states in Table 
EH In Ref. @, .Dj(2580) is identified as 0”, £}(2650) 
is identified as 1~, but one can see that Dj (2580) is 
not fully compatible with the assignment 0“ in our cal¬ 
culation. In this work, we take a different assignment, 
i.e., £j(2580) is identified as dominantly |2 3 S'i) state 
with J p = 1 _ , £}(2650) is identified as |2 3 P 0 ) state 
with J p = 0 + , since our prediction for |2 1 S'o) and |2 3 Si) 
state are lower than the predictions of Godfrey and Isgur 
[lj]. The £j(3000) resonance is identified as dominantly 
|2 3 P2) state with J p = 2 + in this work. 

Since 2004 many new D s mesons have been observed, 
such as £> s j(2632) EH. D* sJ ( 2860) 1, D* sl ( 2710) @, 
P s j(3040) [J|. One can see that our predicted mass 
spectrum for D s mesons below 2600 MeV accommodates 
the experimental data successfully [16], here we discuss 
the assignment of the recently observed D s meson reso¬ 
nances. The D*j(2860) resonance parameters measured 
in Ref. 0] are 

M(D* sJ { 2860)+) = 2862 ±2t 3 MeV 
T = 48 ± 3 ± 6 MeV 


The J p = 0+ assignment is proposed in Ref. 3lj . but the 
observation of P>*j(2860) + —>• D*K decay rules out the 
J p = 0+ assignment [lj. Our mass prediction for 11 3 ZD 3 ) 
state is 2853 MeV, very close to the mass of P>*j(2860), 
thus we suggest that £* ,(2 860) is J p = 3“ state, which 
is also proposed in Ref. |32j . For the £^(2710) state, it 
is measured to be 2710 ± MeV in Ref. fj, the J p = 


1 assignment is proposed in Refs. 


, |3J] , comparing 
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the experimental mass of D* 1 (2710) with our calculation 
for |2 3 5 i) state in Table HU our prediction favors this 
assignment. The meson D s j{ 2632) was firstly observed 
in Ref. [l[ at a mass of 2632.5 ± 1.7 MeV, in our predicted 
mass spectrum for D s mesons we assign it as the |2 3 Po) 
state with J p = 0”, it is the closest states to 2632 MeV, 
although 50 MeV lower. The D s j(3040) resonance is 
observed in the D*K mass spectrum at a mass of 3044 ± 
8^ 30 in Ref. [3j. The assignment of this resonance is 
discussed in detail in Ref. [35J. We identify it as |2 1 Pi) 
with J p = 1 + , since its predicted mass 3012 MeV is the 
closest to the mass of the P s j(3040) resonance. 

For b-flavored meson sector, one can see the predicted 
mass spectra in Table m and m are in good agreement 
with experimental results. Although experimental in¬ 
formation for b-flavored meson spectrum is limited, the 
measured J p = 0 _ , 1~, 1 + and 2 + states can be an extra 
test to our model besides the charmed meson sector. 

Next we would like to discuss the properties of the 
heavy-light quark-antiquark bound states from the point 
view_of heavy quark symmetry. As discussed in Refs. 
EH El, the spin of the heavy qurk decouples from the 
light degrees of freedom in the heavy quark limit. The 
total angular momentum of the light quark is Si = s q + l, 
where s q is the spin of the light quark, l the orbital an¬ 
gular momentum of the light quark. The total angular 
momentum of the heavy light meson is J = s; + sq, 
where sq is the spin of the heavy quark. Since the prop¬ 
erties of the hadronic states do not depend on the spin 
sq and the flavor of the heavy quark due to the heavy 
quark symmetry, the heavy-light meson states only de¬ 
pend on the quantum numbers Si and l , thus we can 
classify the heavy-light mesons by these two quantum 
numbers. Heavy-light meson states can be divided into 
different doublets according to this classification. A dou¬ 
blet with J p = (0“, 1 — )i /2 exists for the case l = 0, for 
example. Mesons within the same doublet degenerate in 
the heavy quark limit. Since the parity of the meson is 
P = (—l) i+1 , the parity P can denote the orbital an¬ 
gular momentum l of the meson in a doublet. For each 
meson state with specific quantum numbers classified in 
the heavy quark limit, we can denote them by | J p )- 

By analyzing angular momentum addition with the 
help of Clebsch-Gordan coefficients, we can decompose 
the basis states | J p ) in the heavy quark symmetry into 
combination of states with definite quantum numbers l , 
S and J, i.e. the state \ 2 S+ 1 Lj). The decompositions of 
all the doublets involved in this work are as follows: 

(1) For l = 0 states in the doublet J p = (0 _ , l~)i/ 2 , 
we get 


and J p = (l + , 2 +) 3 / 2 , we get 

|0+ /2 > = | 3 Po>, (51) 

Itf/a) = V^ |1Pl) + Vl |3Pl) ’ (52) 

l^/ 2 > = -Vf |1Pl) + \/I |3Pl) ’ (53) 

| 2 1 /2 ) = I 3 p 2 >- (54) 

(3) For l = 2 states in the doublets J p = (1 — , 2 — ) 3 / 2 
and J p = ( 2 -, 3 -) b/ 2 , we get 

|l-/ 2 ) = | 3 Pi>, (55) 

l 2 3 - /2> = yfr^} + x/fl 3 ^), (56) 

l 2 5 / 2 > = -'/f| 1 ^)+y|| 3 l? 2 ), (57) 

|3 5 -/2> = | 3 A3>. (58) 

(4) For l = 3 states in the doublets J p = (2 + ,3 + ) 5 / 2 
and J p = (3 + ,4+) 7/2 , we get 

|2 1 /2 ) = 1 3 P 2 ), (59) 

l 3 ^/ 2 > = y|rf 3 )+ 2 ^| 3 P 3 ), (60) 

l 3 7 / 2 > = + y|| 3 P 3 ), (61) 

|4+ 2 ) = | 3 P 4 ). (62) 

(5) For l = 4 states in the doublets J p = (3 - ,4 “) 7 / 2 
and J p = (4“,5“) 9/2 , we get 

l 3 V 2 > = l 3 G 3 ), (63) 

l 4 7 / 2 ) = + -^-| 3 G4), (64) 

l 4 g 72 ) = + gl 3 ^ 4 ), (65) 

|5 9 - /2 > = | 3 G 5 ). (66) 


Particularly, for P = (—1 ) J+1 unnatural parity states , 
they satisfy the relations 

l ^- 1 /2 ) = \/ST |1Zi) + \ZI+T |3 ^’ (67) 

l^/ 2 > = + ( 08 ) 


|0r /2 ) = 1%), (49) 

|ir /2 ) = i 3 ^)- (50) 

(2) For l = 1 states in the doublets J p = (0 + , l + )i /2 


Consequently, the physical states with unnatural par¬ 
ity are linear combinations of I 1 !;) and | 3 ^z), the mixing 
can be described by a mixing angle. We must note here 
that the si gn o f the mixing angle is fraught with ambi¬ 
guities 36, 37|. It depends on the coupling order of the 
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angular momentum and spins, l 7 s\ and s* 2 - The sign 
of the mixing angle will flip if the coupling order of the 
quark spins changes, for example. 

In this work, the properties of the heavy flavored 
mesons implied by the heavy quark limit are approx¬ 
imately maintained in the calculation in the relativis¬ 
tic potential model. The mixing angle in the heavy 
quark limit should be — ArcSin-^/Z/2i + 1. For l = 1, 
we get —ArcSin^/1/3 = —0.615 rad, for l = 2 we have 
—ArcSiny / 2/5 = —0.685 rad, etc. The mixing angle 
increases as l increases. Comparing with the column 
“Multiplet” in the Table I^IV for each meson, one can 
find that this property is maintained in the numerical 
results. As discussed in our previous work []1|, only 
the mixing between 11 1 AA) and |l 3 Pi) in 1 + state of 
cq is very small, which seriously deviates from the re¬ 
sult in the heavy quark limit —0.615 rad, the reason 
has been explained in our previous work [THj]. The 
Belle collaboration determined this mixing angle,which 
is 6 = -0.10 ± 0.03 ± 0.02 ± 0.02 rad Q- Our predic¬ 
tion with small mixing angle for 1 + state of cq is well 


consistent with Belle’s measurement. 

Finally the wave function of each bound state can be 
obtained simultaneously when solving the wave equation. 
Figure Q] is the radial wave functions for D meson as an 
example. It is easy to get all the wave functions when 
they are needed. 


IV Summary 

The bound states of heavy-light quark and antiquark 
system are studied in the relativistic potential model. 
The predictions of heavy-light charmed and 6-flavored 
meson states are presented. The low-lying meson states 
are in good agreement with experiments. The higher 
excitation states are particularly focused on. All the ob¬ 
served mesons can be successfully accommodated in our 
calculation. The wave functions of each bound state can 
be obtained by solving the wave equation and be used in 
the study of B and D decays. 


Acknowledgments 

This work is supported in part by the National Nat¬ 
ural Science Foundation of China under contracts Nos. 
11375088, 10975077, 10735080, 11125525. 


* Electronic address: liujingbin077@mail.nankai.edu.cn 
' Electronic address: yangmz@nankai.edu.cn 

[1] A.V. Evdokimov et al. (SELEX Collaboration), Phys. 
Rev. Lett. 93, 242001 (2004). 

[2] B. Aubert et al. (BaBar collaboration), Phys. Rev. Lett. 
97, 222001 (2006). 

[3] J. Brodzicka et al. (Belle Collaboration), Phys. Rev. Lett. 
100, 092001 (2008). 

[4] B. Aubert et al. (BaBar Collaboration), Phys. Rev. D 80, 
092003 (2009). 

[5] P. del Amo Sanchez et al. (BaBar Collaboration), Phys. 
Rev. D 82, 111101 (2010). 

[6] R. Aaij et al. (LHCb collaboration), JHEP 1309, 145 
(2013). 

[7] A. De R.ujula, H. Georgi and S.L. Glashow, Phys. Rev. 
D12, 147 (1975). 

[8] S. Godfrey and N. Isgur, Phys. Rev. D32, 189 (1985). 


[9] P. Cea, P. Colangelo, G. Nardulli, G. Paiano, and G. 
Preparata, Phys. Rev. D26, 1157 (1982). 

[10] P. Cea, G. Nardulli, and G. Paiano, Phys. Rev. D28, 2291 
(1983). 

[11] P. Cea, P. Colangelo, L. Cosmai and G. Nardulli, Phys. 
Lett. B206, 691 (1988). 

[12] P. Colangelo, G. Nardulli, M. Pietroni, Phys. Rev. D43, 
3002 (1991). 

[13] P. Colangelo, F. De Fazio, M. Ladisa, G. Nardulli, P. 
Santorelli, A. Tricarico, Eur. Phys. J. C8, 81 (1999). 

[14] M.Z. Yang, Euro. Phys. J. C72, 1880 (2012). 

[15] J.B. Liu, M.Z. Yang, JHEP 1407, 106 (2014). 

[16] J. Beringer et al. (Particle Data Group), Phys. Rev. D86, 
010001 ( 2012 ). 

[17] M. Di Pierro and E. Eichten, Phys. Rev. D64, 114004 

( 2001 ). 

[18] E. Eichten, K. Gottfried, T. Kinoshita, K.D. Lane, and 
T.-M. Yan, Phys. Rev. D17, 3090 (1978) [Erratum: Phys. 
Rev. D 21, 313 (1980)]. 

[19] E. Eichten, K. Gottfried, T. Kinoshita, K.D. Lane, and 
T.-M. Yan, Phys. Rev. D21, 203 (1980). 

[20] T. Barnes, S. Godfrey, and E.S. Swanson, Phys. Rev. 
D72, 054026 (2005). 

[21] N. Isgur and M. B. Wise, Phys. Lett. B232, 113 (1989). 

[22] N. Isgur and M. B. Wise, Phys. Lett. B237, 527 (1990). 

[23] E. Eichten and B. Hill, Phys. Lett. B234, 511 (1990). 

[24] E. Eichten and B. Hill, Phys. Lett. B243, 427 (1990). 







8 


TABLE I: Theoretical spectrum of ( cq ) and ( cs ) bound states mainly with the radial quantum number n = 1. The results in 
the column labeled “GI” are theoretical masses from Ref. ;.8j. 



Meson 


Multiplet 

Mass (GeV) 

GI (GeV) 

Exp. (MeV) [6, 16J 


D 

0“ 

MM) 

1.867 

1.88 

1869.62 + 0.15 


D* 

1“ 

0.997|l J Si) - 0.035|1 J A) + 0.051|2Ai) + 0.033I2 13 A) 

—0.026|3 3 S'i> — 0.031|3 3 A) 

2.004 

2.04 

2010.28 + 0.13 




0.046|l J Si) + 0.979|1 J A) - 0.090|2Ai) - 0.148|2' 3 Di) 

—0.041|3 3 S’i) + 0.093|3 3 A) 

2.785 

2.82 



Do(2400) u 

0+ 

fl 3 ^ 

2.302 

2.40 

2318 ± 29 



1+ 

0.017|l i Pi) + 0.990|l 3 Pi> - 0.076|2 i Pi) - 0.117|2 J Pi) 

2.415 

2.44 



L>i(2420) 


0.996|l 1 Pi) - 0.028|l d Pi) - 0.028|2 i Pi) - 0.076|2' 5 Pi) 

2.429 

2.49 

2421.3 +0.6 


D 2 (2460) 

2+ 

0.988|l d P 2 ) - 0.014|l a P 2 ) + 0.135|2 J P 2 ) + 0.011|2 J P 2 ) 

—0.071|3 3 P 2 ) - 0.009|3 3 P 2 ) 

2.468 

2.50 

2464.4 ± 1.9 




0.017|l d P 2 ) + 0.997|l a P 2 ) - 0.024|2 J P 2 ) - 0.047|2 J F 2 ) 

—0.021|3 3 P 2 ) + 0.039|3 3 P 2 ) 

3.123 




Dj(2740)° 

2" 

076911^0^ - 0.639|1 J A) + 0.0003|2 1 P 2 ) + 0.007|2 J A) 

-0.007|3 1 P 2 ) + 0.002|3 3 A) 

2.737 


2737.0 + 3.5 + 11.2 




0.639I1 1 A) + 0.768|1 J A) - 0.005I2 1 A) - 0.021|2 3 A) 

0.013|3 1 P 2 ) + 0.022 3 3 P 2 ) 

2.834 



(cq) 

D}(2760)° 

3" 

0.997|G A) - 0.008|1 J G 3 ) + 0.063|2 J A) + 0.006|2 a G3) 
-0.039|3 3 A) + 0.005|3 3 G 3 > 

2.754 


2760.1 ± 1.1 ±3.7 




0.008|1 J A) + 0.9997|1 j G 3 ) - 0.007|2 J A) + 0.002|2 J G 3 ) 
-0.009|3 3 A) + 0.018|3 3 G 3 ) 

3.397 




D}( 3000)° 

3+ 

0.759|1 1 P 3 > - 0.650|1A 3 ) - 0.016|2 A p3> + 0.017|2 J P 3 ) 

—0.001|3 1 P 3 ) -0.001|3 3 P 3 ) 

2.997 


3008.1 +4.0 




0.65011^3} + 0.759|GP 3 ) + 0.020|2 1 P 3 ) + 0.016|2 J P 3 ) 

+0.002|3 1 P 3 ) + 0.006|3 3 P 3 ) 

3.138 





4 + 

0.9995|1 J P 4 ) - 0.005IGP4) + 0.021|2 J P 4 ) + 0.004|2 3 P r 4 > 

—0.022|3 3 P 4 ) - 0.003|3 3 P 4 ) 

3.003 






0.005|1 J P 4 ) + 0.9994|1 j P 4 ) - 0.003|2 J P 4 ) + 0.031 |2' 3 P/ 4 > 
-0.007|3 3 P 4 ) + 0.008|3 3 P 4 ) 

3.638 





4' 

0.750|1 1 G 4 ) - 0.660|1 j G 4 ) - 0.028|2 1 G 4 ) + 0.026|2 3 G 4 ) 

+0.003|3 1 G 4 ) - 0.003|3 3 G 4 ) 

3.229 






0.660|1 1 G 4 ) + 0.750|1 j G 4 ) + 0.034|2 1 G 4 ) + 0.032|2 3 G 4 > 

—0.002|3 1 G 4 ) + 0.0002|3 3 G 4 ) 

3.398 




Df 

0" 

MM> 

1.963 

1.98 

1968.49 + 0.32 


D? 

1~ 

0.996|l 3 Si) - 0.045|l d A) + 0.052|2Ai) + 0.041|2 J A) 

—0.028|3 3 iS’ 1 ) - 0.036|3 3 A) 

2.103 

2.13 

2112.3 ±0.5 




0.064|l J Si) + 0.950|l d A) - 0.191|2Ai) - 0.201|2 J A) 

—0.047|3 3 S’i) + 0.117|3 3 A) 

2.803 

2.90 



T> s *o(2317) u 

0+ 

MM) 

2.317 

2.48 

2317.8+0.6 


An (2460) 

1+ 

0.475|l ± P > i> + 0.851|1 j Pl) - 0.117|2 1 Pi> - 0.193|2 d Pi) 

2.436 

2.53 

2459.6 + 0.6 


D sl (2536) 


0.871 lM) - 0.488 l d Pi) + 0.040|2 i Pi) - 0.033 2 a Pi) 

2.533 

2.57 

2535.12 +0.13 


AM2573) 

2+ 

0.983|1 J P 2 ) - 0.018|l a P 2 ) + 0.164|2 j P 2 ) + 0.013|2 3 P 2 ) 

—0.080|3 3 P 2 ) - 0.010|3 3 P 2 ) 

2.573 

2.59 

2571.9 + 0.8 




0.027|1 J P 2 ) + 0.991|GP 2 ) - 0.063|2 J P 2 ) - 0.097|2 3 P 2 ) 

—0.036|3 3 P 2 ) + 0.058|3 3 P 2 ) 

3.145 





2” 

0.740|1 1 P 2 ) - 0.672|1 J A) + 0.027|2 1 P 2 ) - 0.007|2•’A) 

-0.017|3 1 P 2 ) + 0.005|3 3 A) 

2.835 






0.671|1 1 P 2 ) + 0.737|1 j P 2 ) - 0.037|2 1 P 2 ) - 0.063|2 3 A) 

+0.025|3 1 P 2 ) + 0.039|3 3 A) 

2.864 



(cs) 

AM2860) 

3” 

0.994M A) - 0.009MA) + 0.093|2 J A) + 0.007|2 3 G 3 ) 

—0.047|3 3 A) - 0.005|3 3 G 3 ) 

2.853 


2862 ± 2±1 [4] 




0.011MA) + 0.998MG 3 ) - 0.014|2 J A) - 0.043|2' 3 G 3 ) 
-0.010|3 3 A) + 0.031|3 3 G 3 ) 

3.416 





3+ 

0.762|1 1 P 3 ) - 0.648|GP 3 ) + O.OOO^ 1 ^) - 0.002|2 3 P 3 ) 

—0.008|3 1 P 3 ) + 0.005|3 3 P 3 ) 

3.092 






0.64811^3) + 0.761|GP 3 ) - 0.00812^) - 0.016|2 3 P 3 ) 

+0.010|3 1 P 3 ) + 0.015|3 3 P 3 ) 

3.165 





4 + 

0.998|1 3 P 4 ) - 0.006|1 j P 4 ) + 0.051|2 3 P 4 ) + 0.004|2 3 P r 4 ) 

—0.030|3 3 F 4 ) -0.003|3 3 F 4 ) 

3.097 






0.006|1 J P 4 ) + 0.9997|1 j P 4 ) - 0.006|2 J F 4 ) - 0.010MA) 
-0.005|3 3 P 4 ) + 0.018|3 3 F 4 ) 

3.653 





4” 

0.754|1 1 G 4 ) - 0.657|1’ 3 G 4 ) - 0.005|2 1 G 4 ) + 0.006|2 3 G 4 ) 

—0.004|3 1 G 4 ) + 0.002|3 3 G 4 ) 

3.320 






0.657|1 1 G 4 ) + 0.754|1 ■ 5 G 4 ) + 0.007|2 1 G 4 ) + 0.004|2 3 G 4 ) 

+0.004|3 1 G 4 ) + 0.007|3 3 G 4 ) 

3.420 
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TABLE II: Theoretical spectrum of (eg) and (cs) bound states mainly with the radial quantum number n = 2. 




Multiplet 

Mass (GeV) 

GI (GeV) 

Meson 

Exp. (MeV) 


0“ 

[2 T Sb5 

2.483 

2.58 




1“ 

—0.043|l J S'i) + 0.083|l J Di) + O^l^i) - 0.055|2 a Di) 
+0.065|3 3 Si) + 0.067|3 3 Di) 

2.613 

2.64 

Pj(2580)° 

2579.5 ± 3.4 ± 5.5[6] 



—0.040|l 3 S'i) + 0.161|l 3 Di) + O.OSS^^Si) + Dx) 

—0.180|3 3 Si> -0.180|3 3 Di) 

3.155 





0+ 

|2°-Po) 

2.685 


P>}(2650) u 

2649.2 + 3.5 + 3.5 [6] 


1+ 

0.083|1 1 P l ) + 0.137|l 3 Pi) + 0.475512^) + 0.867|2 d Pi) 

2.865 






—0.010|l 1 Pi) + 0.014|l J Pi) + 0.887|2 1 Pi) - 0.462|2 d Pi) 

2.902 





2+ 

-0.111|1 3 P 2 ) + 0.029|1 3 P 2 ) + 0.959|2 d P 2 ) - 0.021|2 d P 2 ) 
+0.258 3 3 P 2 ) + 0.016|3 3 P 2 ) 

2.969 


Pj(3000)° 

2971.8 + 8.7 [6] 



-0.057|1 3 P 2 ) + 0.040 1 3 P 2 ) + 0.124|2 d P 2 ) + 0.898|2 d P 2 ) 
-0.412 3 3 P 2 > - 0.059|3 3 P 2 ) 

3.458 




(eg) 

2~ 

—0.003|1 1 P 2 ) - 0.011|1 3 P 2 ) + 0.757|2 i P 2 ) - 0.653|2 3 P 2 ) 
+0.008|3 1 P 2 > + 0.005|3 3 P> 2 ) 

3.144 






0.009|1 1 P 2 ) + 0.020|1 3 P> 2 ) + 0.653|2 i P 2 ) + 0.755|2 d P> 2 ) 
-0.021|3 1 P 2 ) - 0.044|3 3 P> 2 ) 

3.211 





3~ 

-0.057|l d P> 3 ) + 0.009|1 3 G 3 ) + 0.990|2 3 P> 3 ) - 0.010|2 3 G 3 ) 
+0.124|3 3 P> 3 ) + 0.008 3 3 G 3 ) 

3.175 






-0.007|l d P> 3 ) - 0.002|1 3 G 3 ) + 0.012|2 3 P> 3 ) + 0.9996|2 3 G 3 ) 

—0.018|3 3 P> 3 ) -0.015|3 3 G 3 ) 

3.704 





3+ 

0.145|l i P 3 ) - 0.020|1 3 P 3 ) + 0.755|2 i P 3 ) - 0.654|2 3 P 3 ) 

—0.017|3 1 P 3 ) + 0.020|3 3 P 3 ) 

3.360 






—0.018|l i P 3 ) - 0.018 1 3 P 3 ) + 0.654|2 i P 3 ) + 0.755|2 3 P 3 ) 
+0.020|3 1 P 3 ) + 0.011|3 3 P 3 ) 

3.480 





4+ 

-0.020|1 3 P 4 ) + 0.004|1 3 P 4 ) + 0.998|2 3 F 4 ) - 0.006|2 3 P 4 > 
+0.050|3 3 P 4 ) + 0.005|3 3 F 4 ) 

3.373 






—0.004|1 3 P 4 ) - 0.032|1 3 P 4 ) + 0.006|2 3 P 4 ) + 0.999|2 3 P 4 ) 
-0.006|3 3 P 4 ) + 0.029|3 3 F 4 ) 

3.926 





0” 

[Fsb) 

2.582 

2.67 

D sJ ( 2632) 

2632.5 ± 1.7 [1] 


1“ 

-0.035|l d Si) + 0.175|l J Pi) + 0.976|2 d Si) - 0.087|2 d Di) 
+0.063|3 3 Si) + 0.065|3 3 P>i) 

2.709 

2.73 

-Dai (2710) 

2708 ± g^ 11 [3] 



-0.050|l J Si) + 0.190|l J Pi) + 0.082|2 d Si) + 0.816|2 d PL) 
-0.507|3 3 Si> -0.180|3 3 Pi) 

3.191 





0+ 

WPo) 

2.707 





1+ 

0.117|l 1 Pi) + 0.193|l J Pi) + 0.498|2 i Pi) + 0.837|2 d Pi) 

2.911 






—0.044|l 1 Pi) + 0.026|l J Pi) + 0.858|2 1 Pi) - 0.510|2 d Pi) 

3.012 


PL j (3040) 

3044±8+5 U [4] 


2+ 

-0.130|1 3 P 2 ) + 0.070|1 3 P 2 ) + 0.941|2 3 P 2 ) - 0.030|2 3 P 2 ) 

+0.301 3 3 P 2 ) + 0.022|3 3 P 2 ) 

3.070 






-0.001|1 3 P 2 ) + 0.108 1 3 P 2 > - 0.011|2 3 P 2 ) + 0.976|2 3 P 2 ) 

+0.119 3 s P 2 ) -0.143|3 3 P 2 ) 

3.477 




(cs) 

2“ 

0.048|1 1 P 2 ) + 0.041|1 3 P 2 ) - 0.087|2 i P 2 ) + 0.990|2 3 P> 2 ) 

-0.072|3 1 P 2 ) - 0.051|3 3 P> 2 ) 

3.240 






0.011|1 1 P 2 ) + 0.051|1 3 P 2 ) + 0.992|2 i P 2 ) + 0.080|2 3 P> 2 ) 

-0.009|3 1 P 2 ) - 0.080|3 3 P> 2 ) 

3.246 





3“ 

-0.083|1 3 P 3 ) + 0.015|1 3 G 3 ) + 0.982|2 3 P 3 ) - 0.013|2 3 G 3 ) 

+0.171|3 3 P 3 ) + 0.010 3 3 G 3 ) 

3.272 






-0.011|1 3 D 3 ) + 0.044|1 3 G 3 ) + 0.021|2 3 P 3 ) + 0.995|2 3 G 3 ) 

—0.051|3 3 P 3 ) - 0.071|3 3 G 3 ) 

3.724 





3+ 

—0.007|1 1 P 3 ) + 0.0002|1 J P 3 ) + 0 . 754 | 2 1 P r 3 ^ > o.657|2 d F 3 ) 

+0.016|3 1 P 3 > - 0.007|3 3 P 3 ) 

3.454 






0.010|l i P 3 ) + 0.016|l d F 3 ) + 0.657|2 i P 3 ) + 0.753|2 d P 3 ) 

—0.019|3 1 P 3 ) — 0.033|3 3 P 3 ) 

3.507 





4+ 

—0.048|1 J P 4 ) + 0.006|l J iP 4 > + 0.994|2 d P 4 ) - 0.007|2 a P 4 ) 
+0.096|3 3 P 4 ) +0.006|3 3 P 4 ) 

3.467 






-0.005|l d F 4 ) + 0.011|1' 3 P 4 ) + 0.008|2 d P 4 ) + 0.9994|2 J P 4 ) 
-0.012|3 3 P 4 ) -0.027|3 3 F 4 ) 

3.941 












































10 


TABLE III: Theoretical spectrum of ( bq ) and (bs) bound states mainly with the radial quantum number n = 1. 



Meson 


Multiplet 

Mass (GeV) 

GI (GeV) 

Exp. (MeV) [16J 


B 

0~ 

II 1 *} 

5.271 

5.31 

5279.25 +0.17 


B* 

1" 

0.9996|l 3 Si) - 0.012|l’ 3 Di) + O.OlO^Si) + 0.012|2' 3 Di) 
-0.009|3 3 Si> -0.011|3 3 Di) 

5.317 

5.32 

5325.2 +0.4 




0.014|l J Si) + 0.995|l J Z)i> - 0.020|2' 3 S , i> - 0.080|2 J Di> 
-0.011|3 3 Si> + 0.051|3 3 Di> 

6.087 





0+ 


5.696 




Bi(5721) 

1+ 

0.51111^1) + 0.847|l 3 Pi) - 0.082|2 i Pi) - 0.117|2 a Pi) 

5.732 


5723.5 + 2.0 




0.855 ^Pl) - 0.517 l 3 Pi) + 0.031 2 i Pi> - 0.031|2 a Pi) 

5.759 




B 2 *(5747) 

2+ 

0.995|1 3 P 2 ) - 0.005|1 j F 2 ) + 0.090|2 J P 2 ) - 0.00412^2) 
-0.049|3 3 P 2 ) - 0.004|3 3 P 2 ) 

5.775 

5.8 

5743 ± 5 




0.006|1 3 P 2 ) + 0.9995|1 3 P 2 ) - 0.008|2 a P 2 ) - 0.02112^+2) 
-0.006|3 3 P 2 ) + 0.021|3 3 P 2 ) 

6.374 





2” 

0.76411^^2) - 0.645|1 3 P> 2 ) + O.OlO^Da) - 0.005|2 a P 2 > 

-0.010|3 1 P 2 ) + 0.007|3 3 9 2 } 

6.051 






0.644|l i P 2 ) + 0.763|1 3 P> 2 ) - 0.017|2 i P 2 > - 0.025|2 a P 2 ) 

0.016|3 1 P 2 ) +0.021 3 3 P 2 ) 

6.105 



( 6 <?) 


3” 

0.999|1 3 D 3 ) - 0.003|1 3 G 3 ) + 0.038 2 J P> 3 ) + 0.002|2 J G 3 ) 
-0.025|3 3 P> 3 ) - 0.002|3 3 G 3 ) 

6.061 






0.003|1 3 P 3 ) + 0.9999|1 J G 3 > - 0.003|2 a P 3 ) + 0.007|2 J G 3 ) 
-0.002|3 3 P> 3 ) + 0.010|3 3 G 3 ) 

6.614 





3+ 

0.753|l i P 3 ) - 0.658|1 J P 3 ) - 0.006|2 3 P 3 ) + OW^Pa) 

—0.004|3 1 P 3 ) + 0.003|3 3 P 3 ) 

6.296 






0.658|l i P 3 ) + 0.753|1 3 P 3 ) + OW^Ps) + O.OOS^Pb) 

+0.005|3 1 P 3 ) + 0.007|3 3 P 3 ) 

6.383 





4+ 

0.9999|1 3 P 4 ) - 0.002|l l3 P 4 > + 0.011|2^P 4 ) + 0.001|2 a P 4 ) 

—0.013|3 3 P 4 ) - 0.001|3 3 P 4 ) 

6.302 






0.002|1 J P 4 ) + 0.9997|1 J P 4 ) - 0.002|2^P 4 ) + 0.024|2 a P 4 ) 

—0.001|3 3 P 4 ) + 0.004|3 3 P 4 ) 

6.826 





4” 

0.744|1 1 G 4 ) - 0.667|1 3 G 4 ) - 0.017|2 1 G 4 ) + 0.016|2 ;> G 4 ) 

—0.0004|3 1 G 4 ) -0.0001|3 3 G 4 ) 

6.511 






0.667|1 1 G 4 ) + 0.744|1 J G 4 ) + 0.020|2 1 G 4 ) + 0.020|2 ;> G 4 ) 

+0.001|3 1 G 4 ) +0.002|3 3 G 4 ) 

6.619 




B s 

0“ 

II 1 *} 

5.360 

5.39 

5366.77 ± 0.24 


B* s 

1" 

0.9995|l J S'i) - 0.014|l' 3 Di) + 0.020|2' 3 Si) + O.OlS^Pi) 
-0.010|3 3 Si) -0.012|3 3 Di) 

5.408 

5.45 

5415.4± 3 ;! 




0.017|l J Si) + 0.989|l 3 P>i) - O.OSO^Si) - 0.122|2 J Di) 
-0.010|3 3 Si) - 0.069|3 3 Pi) 

6.120 





0+ 


5.731 





1+ 

0.53711^1) + 0.822|l J Pi) - 0.108|2 i Pi) - 0.159|2 a Pi) 

5.765 




B s i(5830) 


0.835 ^Pl) - 0.546 l 3 Pi) + 0.056 2 i Pi) - 0.042|2 a Pi) 

5.850 


5829.4 + 0.7 


B s * 2 (5840) 

2+ 

0.992|1 J P2> - 0.006|1 3 P 2 ) + 0.112|2 J P 2 ) + O.OOS^Pa) 
-0.056|3 3 P 2 ) - 0.004|3 3 P 2 ) 

5.866 

5.88 

5839.7 + 0.6 




0.008|1 3 P 2 ) + 0.998|1 J P 2 ) - 0.015|2 J P 2 ) - O^^Pa) 

—0.008|3 3 P 2 ) + 0.033|3 3 P 2 ) 

6.405 





2~ 

-0.684|l i P 2 > + 0.728|1 3 P> 2 ) - 0.033|2 i P> 2 ) + 0.011|2 J P> 2 ) 

+0.019 3 ^ 2 ) - 0.007 3 3 P> 2 ) 

6.134 






0.727|l i P 2 ) + 0.682|l a P> 2 ) - 0.037|2 i P 2 > - 0.058|2 a P 2 ) 

+0.023|3 1 P 2 ) + 0.034|3 3 9 2 } 

6.142 



(bs) 


3” 

0.998|1 3 D 3 ) - 0.004|PG 3 ) + 0.061|2 J P> 3 ) + 0.003|2 J G 3 ) 
-0.031|3 3 P> 3 ) - 0.002|3 3 G 3 ) 

6.144 






0.004|1 J P> 3 ) + 0.9995|1 3 G 3 ) - 0.005|2 J P> 3 ) - 0.024|2 J G 3 ) 
-0.003|3 3 P> 3 ) + 0.018|3 3 G 3 ) 

6.641 





3+ 

0.749|l i P 3 ) - 0.662|1 j P 3 ) + O.Oll^Pj) - 0.007|2 3 F 3 ) 

—0.009|3 1 P 3 ) + 0.007|3 3 P 3 ) 

6.374 






0.662|l i P 3 ) + 0.749|1 j P 3 ) - O.OM^Pa) - O.OIO^Pb) 

+0.011|3 1 P 3 ) + 0.014|3 3 P 3 ) 

6.416 





4+ 

0.9992|1 3 P 4 ) - 0.002|1 3 P 4 ) + 0.033|2 3 P 4 ) + 0.002|2 3 P 4 ) 

—0.019|3 3 P 4 ) - 0.001|3 3 P 4 ) 

6.380 






0.002|1 3 P 4 ) +0.9999|1 3 P 4 ) - 0.002|2 3 P 4 ) - 0.006|2 3 P 4 ) 

-0.001|3 3 P 4 ) + 0.010|3 3 F 4 ) 

6.849 





4” 

0.744|1 1 G 4 ) - 0.668|1 3 G 4 ) + 0.001|2 1 G 4 ) + 0.001|2 3 G 4 ) 

—0.005|3 1 G 4 ) + 0.004|3 3 G 4 ) 

6.585 






0.668|1 1 G 4 ) + 0.744|1 3 G 4 ) - 0.001|2 1 G 4 ) - 0.003|2 3 G 4 ) 

+0.006|3 1 G 4 ) +0.007|3 3 G 4 ) 

6.647 




























































11 


TABLE IV: Theoretical spectrum of ( bq ) and (6s) bound states mainly with the radial quantum number n = 2. 




Multiplet 

Mass (GeV) 

GI (GeV) 


0” 

r I^So) 

5.836 

5.90 


1” 

-0.018|T 3 Si) + 0.019|V_Di) + 0.999|2' 3 S'i) - 0.016|2 a P>i) 
+0.022|3 3 5 'i) + 0.014|3 3 Di) 

5.878 

5.93 



-0.034|l l5 S'i) + 0.086|l 3 Di) + 0.016|2' 3 S'i) + 0.988|2 a L>i) 
-0.034|3 3 S , i) -0.120|3 3 Di) 

6.422 



0+ 

|2 J P 0 > 

6.070 



1+ 

r O.OSlIEPi) + 0.118|l a Pi) + 0 . 55212 ^ 1 } + 0.822|2 J Pi) 

6.137 




—0.036|l 1 Pi) + 0.025| VPi) + 0.829|2 i Pi) - 0.557|2 J Pi) " 

6.196 



2+ 

-0.080|l a P 2 } + 0.009|1 J P 2 ) + 0.98312^2) - 0.007|2' 3 E 2 ) 
+0.168|3 3 P 2 ) + 0.006 3 3 P 2 ) 

6.220 




-0.007|l a P 2 } + 0.022|1 j P 2 ) + 0.012|2 ^ P 2 ) + 0.998|2’ 5 P 2 ) 

—0.031|3 3 P 2 ) - 0.044 3 3 P 2 ) 

6.676 


(bq) 

2~ 

-0.011|l i P 2 ) + 0.004|1 J D 2 ) + 0.753 2 L D 2 ) - 0.657|2 J D 2 ) 
+0.024 \3 1 D 2 ) - 0.014 3 3 P> 2 ) 

6.417 




0.018|l i P 2 ) + 0.026|1 J D 2 ) + 0.656|2 i P 2 ) + 0.752|2 d D 2 ) 
-0.034|3 1 P 2 ) - 0.048|3 3 P> 2 ) 

6.446 



3" 

-0.037|1 j P> 3 ) + 0.003|1 J G 3 ) + 0.996|2 3 .D 3 ) - 0.004|2 J G 3 ) 
+0.076|3 3 P> 3 ) + 0.003 3 3 G 3 ) 

6.430 




-0.003|1 J P> 3 ) - 0.007|1 J G 3 ) + 0.004|2 J P> 3 ) + 0.9999|2 3 G 3 > 
-0.006|3 3 P> 3 ) - 0.002|3 3 G 3 ) 

6.891 



3+ 

0.00611^3) - 0.007|1 J P 3 ) + 0.75112^3) - 0.660|2 a P 3 ) 

— 0.002|3 1 P 3 ) + 0.004|3 3 P 3 ) 

6.618 




— 0.007|l i P 3 ) - 0.006|1 j P 3 ) + 0.660 2 ^ 3 ) + 0.751|2 3 P 3 ) 
+0.002|3 1 P 3 ) - 0.002 3 3 P 3 ) 

6.689 



4+ 

-0.010|1 J P 4 ) + 0.002|l a P 4 ) + 0.9995|2 J P 4 ) - 0.002|2 a P 4 ) 
+0.028|3 3 P 4 ) + 0.002|3 3 P 4 ) 

6.626 




-0.002|1 J P 4 ) - 0.024|1 J P 4 ) + 0.002|2 3 F 4 ) + 0.999|2 J P 4 ) 
-0.003|3 3 P 4 ) + 0.025|3 3 P 4 > 

7.085 



0" 

WSo) 

5.928 

5.98 


ir 

-0.018|l 3 5i) + 0.028|l 3 Pi) + O.ggg^ 3 ^!) - 0.019|2 3 P>i) 
+0.023|3 3 Si) + 0.016|3 3 Pi) 

5.970 

6.01 



-0.014|l 3 5i) + 0.132|l 3 P>i) + O.OlS^Si) + 0.976|2 3 P>i) 
-O.OSSIS 3 *) -0.160|3 3 Pi) 

6.460 



0+ 

|2 J P 0 ) 

6.109 



1+ 

O.lOSlPPi) + 0.159|l 3 Pi) + 0.547|2 ± P > i> + 0.815|2 3 Pi) 

6.187 




—0.058|l 1 Pi) + 0.039|l a Pi) + 0.828|2 1 Pi) — 0.556|2 J Pi) 

6.289 



2+ 

—0.098|1' 3 P 2 ) + 0.016|1 J P 2 ) + 0.974|2 15 P 2 ) - 0.009|2 3 F 2 ) 
+0.202|3 3 P 2 ) + 0.007 3 3 P 2 ) 

6.308 




—0.024|1' 3 P 2 ) + 0.056|1 J P 2 ) + 0.055|2' 3 P 2 ) + 0.965|2 3 F 2 ) 

—0.235|3 3 P 2 ) - 0.081 3 3 P 2 ) 

6.709 


(1bs ) 

2" 

0.046|l i P 2 ) + 0.057|1 J P 2 ) + 0.561|2 i P 2 ) + 0.818|2 d P> 2 ) 
-0.066|3 1 P 2 ) - 0.080|3 3 P> 2 ) 

6.483 




-0.024|l i P 2 ) + 0.022 1 J Z) 2 ) + 0.822 2 L D 2 ) - 0.565|2 J D 2 ) 
+0.044 3 1 P 2 ) - 0.042 3 3 P> 2 ) 

6.499 



3" 

-0.057|1 J P> 3 ) + O.OOSll^Gs) + 0.992|2 a D 3 ) - 0.005|2 J G 3 ) 
+0.109|3 3 P> 3 ) + 0.004 3 3 G 3 ) 

6.512 




-0.003|1 J P> 3 ) + 0.025|1' 3 G 3 ) + 0.006|2 3 .D 3 ) + 0.999|2 J G 3 ) 

—0.012|3 3 P 3 ) - 0.044|3 3 G 3 ) 

6.919 



3+ 

—0.011|l i P 3 ) + 0.007|1 j P 3 ) + 0.740|2 i P 3 ) - 0.672|2 J P 3 ) 
+0.023|3 1 P 3 ) - 0.016|3 3 P 3 ) 

6.695 




0.015|l i P 3 ) + 0.020|1 j P 3 ) + 0.651|2 i P 3 ) + 0.740|2 13 P 3 ) 

—0.026|3 1 P 3 ) - 0.035|3 3 P 3 ) 

6.722 



4+ 

-0.032|1 J P 4 ) + 0.003|1 J P 4 ) + 0.998|2' 3 P 4 ) - 0.003|2 J P 4 ) 
+0.061|3 3 P 4 ) + 0.002|3 3 P 4 ) 

6.704 




-0.002|1 j P 4 ) + 0.006|1 j P 4 ) + 0.003|2 J P 4 ) + 0.9998|2 J P 4 ) 
-0.005|3 3 P 4 ) -0.016|3 3 P 4 ) 

7.108 
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